Contact values of the particle-particle and wall-particle correlation functions in a 

hard-sphere polydisperse fluid 

Andres Santofl and Santos B. Yust^ 
Departamento de Fisica, Universidad de Extremadura, Badajoz, E-06071, Spain 

Mariano Lopez de Harc[j] 
Centro de Investigation en Energia, UNAM, Temixco, Morelos 62580, Mexico 
(Dated: February 2, 2008) 

The contact values g(a, a') of the radial distribution functions of a fluid of (additive) hard spheres 
with a given size distribution f(o~) are considered. A "universality" assumption is introduced, 
according to which, at a given packing fraction rj, g(a,a') — G(z(a, a')), where G is a common 
function independent of the number of components (either finite or infinite) and z(a, a') = [2a a 1 / (a+ 
o~')\n2/ M3 i s a dimensionless parameter, [i n being the n-th moment of the diameter distribution. A 
cubic form proposal for the z-dependence of G is made and known exact consistency conditions 
for the point particle and equal size limits, as well as between two different routes to compute the 
pressure of the system in the presence of a hard wall, are used to express G(z) in terms of the 
radial distribution at contact of the one-component system. For polydisperse systems we compare 
the contact values of the wall-particle correlation function and the compressibility factor with those 
obtained from recent Monte Carlo simulations. 



I. INTRODUCTION 

The prominence of hard-core fluids in liquid state the- 
ory as prototype systems for theoretical understanding 
and stepping stone for the study of more realistic fluids 
can hardly be overemphasized. It is well known that the 
form of the pressure equation of these systems acquires a 
particularly simple representation in terms of the contact 
values of the radial distribution functions (rdf). There- 
fore, since for hard-core fluids the internal energy reduces 
to that of ideal gases, the knowledge of such contact val- 
ues would be enough to obtain their equation of state 
(EOS) and all their thermodynamic properties. 1 Unfor- 
tunately up to the present day, and except for the case of 
hard rods, no exact expressions either for the contact val- 
ues of the rdf or for the EOS are available and so, in order 
to make progress, research has relied on various approx- 
imate (mostly empirical or semiempirical) approaches or 
on the results of simulations. The situation is rather more 
complicated for mixtures than for single component flu- 
ids and hence it is not surprising that studies for the 
former are fewer. From the analytical point of view, per- 
haps the most important result is the exact solution of 
the Percus-Yevick (PY) equation of additive hard-sphere 
mixtures carried out by Lebowitz^ that includes explicit 
expressions for the contact values of the rdf. In turn, 
such expressions served as the basis for the derivation of 
the widely used and rather accurate Boublfk-Mansoori- 
Carnahan-Starling-Leland (BMCSL) EOS 3 ' 4 for hard- 
sphere mixtures. In fact, Boubh'ls 3 - (and, independently, 
Grundke and Henderson^ and Lee and Levesqua^) pro- 
posed an interpolation between the PY contact values 
and the ones of the Scaled Particle Theory (SPT)A£ We 
will refer to this interpolation as the Boublik-Grundkc- 
Henderson-Lee-Levesque (BGHLL) approximation for 
the contact values, refinements of which have been subse- 
quently introduced, among others, by Henderson et al.^ 



Matyushov and LadanyijiS and Barrio and Solanaii to 
eliminate some drawbacks of the BMCSL EOS in the so- 
called colloidal limit of binary hard-sphere mixtures. 

It is interesting to point out that in the case of multi- 
component mixtures of hard spheres, the contact values 
which follow from the solution of the PY equation^ those 
of the SPT approximation^^ and those of the BGHLL 
interpolation 3 ^^ present a kind of "universal" behavior 
in the following sense. Once the packing fraction is fixed, 
the expressions for the contact values of the rdf for all 
pairs of like and unlike species depend on the diameters 
of both species and on the size distribution only through 
a single dimensionless parameter, irrespective of the num- 
ber of components in the mixture. In previous workiSi 1 ^ 
we have introduced approximate expressions for the con- 
tact values of the rdf valid for mixtures with an arbitrary 
number of components and in arbitrary dimensionality, 
that require as input the EOS of the one-component fluid. 
Apart from satisfying known consistency conditions, they 
are sufficiently general and flexible to accommodate any 
given EOS for the single fluid and also share the univer- 
sal behavior alluded to above. In the latter paper ji£ two 
functional forms (a quadratic one and a rational one) 
were examined. We found that the best global agree- 
ment with the available simulation results for binary and 
ternary mixtures was provided by the quadratic function, 
which has a structure similar to the SPT and BGHLL 
prescriptions, except perhaps for very disparate mixtures, 
where the rational approximation seemed to be prefer- 
able. 

The universality feature present in the above propos- 
als, which applies to mixtures with an arbitrary number 
of components N and an arbitrary size distribution, per- 
mits in principle to consider different situations. For in- 
stance, one could study the structural properties of an 
M-component mixture in the presence of a hard wall 
by considering a mixture with N = M + 1 components 
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and taking the limit in which the diameter of one of the 
species goes to infinity. Also, one could take the limit 
N — > oo corresponding to a polydisperse system of hard 
spheres in which rather than a discrete set of values for 
the diameters, one has a continuous distribution. Interest 
in studying the thermodynamic and structural properties 
of these polydisperse systems dates back to the late 1970s 
and the 1980si^ and has recently been revived^ Of par- 
ticular concern to us here is a recent paper by Buzzacchi 
et alw& in which they study the structural properties of 
polydisperse hard spheres in the presence of a hard wall. 
It seems natural to compare their results with the ones 
derived from other theories i n particular from 

our approach^Siii However, as shown below, except for 
the SPT contact values (which are known to be gener- 
ally less accurate than other proposals), the rest of the 
approximations (including our previous proposals^^) 
lead to an inconsistency between two different ways of 
computing the pressure in the polydisperse fluid. 

The major aim of this paper is to provide yet another 
(more general) approximation for the contact values of 
the rdf that preserves the property of universality present 
in the PY, SPT, and BGHLL, as well as in our previous 
approximations, but avoids the problem just stated. It 
is with this new approximation that we will compare the 
results of Ref . [l^ 

The paper is organized as follows. In Sec. |n] we de- 
rive the new proposal for the contact values of the rdf 
using the known consistency conditions and two differ- 
ent routes to compute the compressibility factor of the 
polydisperse hard-sphere system. Section IlIII deals with 
the comparison between our contact values, the ensuing 
compressibility factors, and the results of Buzzacchi et 
al*& and other theories. We close the paper in Sec. IIVI 
with further discussion and some concluding remarks. 



II. CONTACT VALUES OF THE RADIAL 
DISTRIBUTION FUNCTIONS 



A. Our proposal 

Let us consider a polydisperse hard-sphere mixture 
with a given size distribution /(cr) (either continuous or 
discrete, the latter being of the form f(a) = J^i x i^( a i ~ 
cr)) at a given packing fraction T) = fp/J,3, where p is the 
(total) number density and 



P 00 

= (a n )= daa n f(a) 
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(1) 



denotes the n-th moment of the size distribution. 

We will use the notation g(a,a') for the contact value 
of the pair correlation function of particles of diameters 
a and & ' . This function enters into the virial expression 



of the EOS asii 

z - P = 
pk B T 



1 + 4- 



da 

P-3 Jo 
J\3 



da'f(a)f(a') 



0(0", cr') 



= l+^{(* + o'yg(v,o')), (2) 

where p is the pressure and T is the absolute temperature. 

A hard wall can be seen as a sphere of infinite diameter. 
As a consequence, the contact value of the correlation 
function g w (cr) of a sphere of diameter a with the wall is 
obtained from g(a,a') as 



g w (a) = lim g(cr,cr'). 



(3) 



An alternative route to the EOS is then provided by us- 
ing the sum rule connecting the pressure and the above 
contact values^ namely 



da f(cr)g w (cr) = (g w (cr)). 



(4) 



The subscript w in Z w has been used to emphasize that 
Eq. represents a route alternative to the virial one, 
Eq. J2J, to get the EOS of the hard-sphere polydisperse 
fluid. Of course, Z = Z w in an exact description, but 
Z and Z w may differ when dealing with approximate ex- 
pressions for g(a,a') and the associated g w {u). 

Now we consider a class of approximations of the 
typi 
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1 (a,a')=G(z(cr,a')), 



where 



z(a, a 



2aa' fi2 
cr + a' n 3 



(5) 



(6) 



is a dimensionless parameter depending on the diame- 
ters a and a', as well as on the second and third mo- 
ments. According to Eq. JSJ, at a given packing fraction 
•q all the dependence of <?(cr, a') on cr, cr', and the details 
of the size distribution /(cr) occurs through the single 
parameter z(a,a'). Once one accepts the "univerality" 
ansatz l|5)l. it remains to propose an explicit form for the 
function G(z). To that end, some consistency conditions 
might be useful. First, in the one-component limit, i.e., 
/(cr) = S(cr — ctq), one has z = 1, so tha1>i2ii^ 



G(z = 1) = g p , 



(7) 



where g p is the contact value of the radial distribution 
function of the one-component fluid at the same packing 
fraction 77 as the packing fraction of the mixture. Next, 
the case of a mixture in which one of the species is made 
of point particles, i.e., a — > 0, leads toi2*H 



G{z = 0) 



1 



1-7? 



Gq. 



(8) 
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Conditions Q and (JSJ are the basic ones. A more strin- 
gent condition is the self-consistency between the routes 
J2J and |@J for any distribution /(er). To proceed further, 
let us express G(z) as a series in powers of z: 



G{z) = G Q + J2 G ^ 



(9) 



where it has been assumed that z = is a regular point. 
Condition © is already built in. In agreement with the 
universality assumption JSJ , the coefficients G n are inde- 
pendent of the size distribution, being functions of the 
packing fraction ry only. After simple algebra, the com- 
pressibility factor obtained by inserting the ansatz JSJ, 
along with Eq. ©, into Eq. (0) reads 



Z = Go 



3?7 G 

A*3 



n=l 



Pi 



x(a a (<j + <t ) ). 



(10) 



At this point, we impose the condition that this com- 
pressibility factor depends functionally on the size dis- 
tribution f(cr) only through a finite number of moments. 
This implies that the series in Eq. (|10|l must be truncated 
after n = 3. Therefore, we restrict ourselves to the class 
of approximations 



G(z) =Go + G lZ + G 2 z 2 + G 3 z A . (11 
Using the approximation (|llfl in Eq. IjlUfl we get 
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Z = Ga+rj (3G + 2G X ) + 2^ (Gi + 2G 2 + 2G 3 ) 

(12) 

Note that the dependence of Z on f(a) through n 2 , 
and /X3 is explicit. It only remains to determine the r\- 
dependence of G\, G2, and G3. 

Now we turn to the alternative route to derive the com- 
pressibility factor using Eq. Q). From Eqs. 10, (JSJ), and 
(|ll|l one obtains the approximation 



g w (a) = G + G 1 z w (a) + G 2 z 2 w {a) + G 3 zl(a), 
where 

Zw{&) = lira z(<j 1 a') = 2a—. 

(j'^OO fj, 3 

Thus the EOS (gj then becomes 
Z w = Gq 



2 ^ Gl+4 4(G 2 
A*3 M3 



2G 3 



(13) 



(14) 



(15) 



Again, the dependence of Z w on the distribution mo- 
ments is explicit. In fact, both Z and Z w are linear in 
the combinations \i\ 1^2/^3 and ^2/^3- The difference be- 
tween Eqs. (|12f) and (|15fl is given by 

ry ry ^1^2 



M3 

-2(1 



[3r?G - 2(1 - ? ? )G 1 ] + 
r,)G 2 - 2(2 - 77)G 3 ] 
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If we want to have Z = Z w for any dispersity, the 
coefficients of ^1^2/^3 and of /J^/mI in Eq. (|16fl must 
vanish simultaneously. This gives 
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3rj 
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G, 



(17) 



(18) 



where we have made use of the definition of Go, Eq. JSJ. 

An extra condition is required to close the problem. 
This follows from the equal size limit given in Eq. J7J, 
which after some algebra yields 



G 2 = (2 - r,)g p 



2 + rf/A 



with 



G 3 = (l-r,) ( 5 f T - ff p) 
jspt _ 1 - »?/2 + ?? 2 /4 



(19) 



(20) 
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the contact value of the radial distribution function of 
the one-component fluid in the SPT. It is also interesting 
to point out that from Eqs. (|19f) and (|2U|I it follows that 



where 



G 2 + G3 



2p PY 



PY 

5 P - 9 P , 



1 + 77/2 

(1 - rj? 



(22) 



(23) 



is the contact value of the rdf for a one-component fluid 
in the PY theory. 

In summary, from Eqs. JSJ), ©, {TTJ, {I7|), {EH, and 
(|20|l we finally get the following expression for the contact 
value of the particle-particle rdf: 
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Similarly, the wall-particle expression is 
1 3?7 /j 2 



(24) 



5iu(o-) 



1 - V (1-r?) A*3 



+4 



(2 - V)9 P 



2 + if /A 



+8(l-7 ? )(. 9 f T - 5p ) ^ 



( 

A*3 

3 



/'3 



(25) 



4 



Approximation 
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TABLE I: Contact values of the one-component fluid g p , coefficients Gi and Gi, and difference Z — Z w corresponding to 
different theories. 



With the above results the compressibility factor may be 
finally written in terms of g p as 
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Equations (|24H are the main results of this paper. 
Note that the contact values of the system, and hence the 
EOS, are wholly determined once the EOS of the one- 
component fluid (and thus g p ) is chosen. Therefore, our 
proposal remains general and flexible in the sense that, 
while fulfilling the consistency conditions 0, JHJ), and 
Z = Z w , the choice of g p can be done at will. Henceforth 
we will denote our approximation by the label "e3" to em- 
phasize that (i) it extends any desired g p to the polydis- 
perse case and (ii) G(z) is a cubic function of z. When a 
particular one-component approximation "A" is chosen, 
we will use the superscript "eA3" to refer to its extension. 
For instance, insertion of the Carnahan-Starling EOS 19 
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into Eqs. and gives g eCSa (a,a') and s° CS3 (cr), 
respectively. More specifically, 
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Mixture 


Type 




fin/ Hi 


M3/V? 


V 


1 


Top-hat (c = 


0.2) 


1.0133 


1.04 


0.2 


2 


Top-hat (c = 


0.2) 


1.0133 


1.04 


0.4 


3 


Top-hat (c = 


0.7) 


1.1633 


1.49 


0.4 


4 


Schulz (q — 


5) 


1.1667 


1.5556 


0.2 


5 


Schulz (q — 


5) 


1.1667 


1.5556 


0.4 



TABLE II: Parameters of the size distributions for the exam- 
ined mixtures. 



B. Connection with former work 

As mentioned in the Introduction, the PY, 2 
the SPT^ the BGHLL, 3 -^ and our previous 
approximations^ 2 *^ for the contact values of the rdf share 
the universality property indicated in Eq. Further- 
more, they have a polynomial dependence on z: linear in 
the case of the PY approximation and the one we pro- 
posed in Ref. (here termed as "el"); quadratic in the 
case of the SPT and BGHLL approximations, as well 
as in our quadratic proposal of Ref. [13| (here termed as 
"e2"). Thus, these five approximations (of which only el 
and e2 allow for a free choice of the one-component rdf 
g p ) may also be expressed in the form of Eq. (|llfl . except 
that G3 = 0. The corresponding coefficients Gi and G2 
appear in Table D Further, we have also included a re- 
cent proposal by Viduna and Smith (VS)S1 that may be 
cast into the form of Eq. Hll|l (again with G3 = 0) but 
does not comply with the ansatz (jSJ) since the coefficient 
Gi depends on the moments of the size distribution. Fi- 
nally we have also included a column with the difference 
between Z and Z w for all those theories. 

From Table [3 we observe that, among the approxima- 
tions with G3 = 0, only the SPT approximation yields 
a consistent EOS through the virial and the hard-wall 
routes, for any density and any degree of polydisper- 
sity. On the other hand, this internal consistency is at 
the expense of the rather poor quality of the SPT con- 
tact value <7p PT in the one-component case. The PY, 
BGHLL, el and VS approximations are not consistent 
with Z — Z w , even in the one-component limit (in which 
case ii n — ► /i"). In the case of the e2 approximation, the 
inconsistency decreases with the degree of dispersity and 
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Mixture 


MC 


z py 


yPY 


z spt 


^BMCSL 


yBGHLL 


^cCSl 


^oCSl 




7-cCS2 


z vs 


yVS 


1 


2.374 


2.344 


2.163 


2.389 


2.374 


2.314 


2.374 


2.240 


2.374 


2.373 


2.376 


2.377 


2 


6.746 


6.197 


4.915 


7.052 


6.767 


6.340 


6.771 


5.636 


6.765 


6.762 


6.792 


6.743 


3 


5.479 


5.215 


4.269 


5.845 


5.635 


5.320 


5.656 


4.848 


5.622 


5.603 


5.663 


5.642 


4 


2.110 


2.076 


1.953 


2.107 


2.097 


2.056 


2.098 


2.012 


2.096 


2.091 


2.099 


2.102 


5 


5.634 


5.042 


4.167 


5.625 


5.431 


5.138 


5.458 


4.722 


5.414 


5.389 


5.461 


5.448 



TABLE III: Comparison between the compressibility factor as obtained by MC simulations 16 and with the different theories. 
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FIG. 1: (Color online). Plot of the difference of contact values 
gw — g w FT as a function of a for a polydisperse mixture with 
a top-hat distribution 13UH with c = 0.2 at a packing fraction 
r/ = 0.2. The symbols are MC simulations.— The lines are 

SPT (■■•), BGHLL ( ), VS ( ), eCS2 ( ), and 

eCS3 (— ). 
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FIG. 2: (Color online). Plot of the difference of contact values 
gw — <?ui PT as a function of a for a polydisperse mixture with 
a top-hat distribution QMIH with c = 0.2 at a packing fraction 
r\ — 0.4. The symbols are MC simulations;— The lines are 

SPT (-••), BGHLL ( ), VS ( ), eCS2 ( ), and 

eCS3 (— ). 



disappears in the one-component limit. 

It must be noted that the el approximation embod- 
ies the PY approximation as a particular case, i.e., the 
choice <7 P = g PY yields G ePY1 = G PY . Analogously, the 
SPT approximation can be recovered from e2 and e3: 
,9 P = gl PT => G cSPT2 = G cSPT3 = G SPT . 

Another comment is in order at this stage. From Eq. 
(|12fl we can observe that, for the class of approxima- 
tions l|llf) . the compressibility factor Z does not depend 
on the individual values of the coefficients G2 and G3, 
but only on their sum. As a consequence, two differ- 



ent approximations of the form (|llfl sharing the same 
density dependence of Gi and G2 + G3 also share the 
same virial EOS. For instance, if one makes the choice 
g p = ,g PY , then Gf Y3 + Gf Y3 = 0, see Eq. (Q), and 
so Z cPY3 = Z PY , even though G cPY3 (z) ^ G PY (z). Fur- 
thermore, if one makes the more sensible choice g p = g^ s , 
then Gf S3 +Gf S3 = Gf GHLL , so that Z cCS3 = Z BMCSL , 
but again G cCS3 (z) ^ G BGHLL (z). 
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FIG. 3: (Color online). Plot of the difference of contact values 
ffiu — <7u; PT as a function of a for a polydisperse mixture with 
a top-hat distribution (1301) with c = 0.7 at a packing fraction 
r/ = 0.4. The symbols are MC simulations.— The lines are 

SPT (■■•), BGHLL ( ), VS ( ), eCS2 ( ), and 

eCS3 (— ). 



FIG. 4: (Color online). Plot of the difference of contact values 
9w — 3™ ?T as a function of a for a polydisperse mixture with 
a Schulz distribution 1311 with q — 5 at a packing fraction 
t] — 0.2. The symbols are MC simulations;— The lines are 

SPT (■■■), BGHLL ( ), VS ( ), eCS2 ( ), and 

eCS3 (— ). 



III. COMPARISON WITH MONTE CARLO 
SIMULATIONS 

Buzzacchi et ali& have recently computed by Monte 
Carlo (MC) simulations the wall-particle contact value 
g w {&) and the compressibility factor Z for polydisperse 
hard spheres of packing fractions r\ — 0.2 and r\ = 0.4 
with either a top-hat distribution of sizes given by 



it what follows, we have made the choice 



/(") = 



l/2c, m (1 - c) < a < /ii(l 



0. 



otherwise 



or with a Schulz distribution of the form 

/(cr) = — a exp a 

qlfii V Mi J V Mi 



(30) 



(31) 



In Table ITT1 we present the values of the parameters cor- 
responding to the examined mixtures. 

Table ITTTl compares the MC data of Z for mixtures 1-5 
with values obtained from different theoretical proposals 
by using both the virial and the wall routes. Here and 



<h 



— ^cs 



the approximations el, e2, and e3. As indicated above, 

Z cCS3 = Z cCS3 = ^BMCSL and ^SPT = ^SPT Notg ^ 
Z BGHLL = ^BMCSL 5 but ^BGHLL ^ ^BMCSL j t can be 

observed from Table II I II that (apart from the SPT and 

the eCS3) the eCS2 and VS expressions for the contact 

values provide the least internal inconsistency between 

Z and Z w . We can also observe that Z BMCSL , Z cCS1 , 

Z oCS2 _ z^ CS2 , and Z vs ~ Z^ s are the most accurate 

EOS for the top-hat mixtures 1-3. On the other hand, 

the most accurate EOS for the Schulz cases 4 and 5 is 
z spt 

Now we turn to the main topic of this paper, namely 
the contact values of the rdf. For the same polydisperse 
systems considered in Table UTI in Figs. ^ El we show the 
comparison between the simulation results^ for the wall- 
particle correlation function g w (a) and those of different 
theories. Since all theories yield contact values which are 
increasing functions of cr, in order to emphasize features 
that would otherwise be difficult to ascertain, we have 
decided to represent the difference g w — g^, PT rather than 
g w . The choice of g£, PT as reference was motivated by the 
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FIG. 5: (Color online). Plot of the difference of contact values 
gw — <7™ PT as a function of o for a polydisperse mixture with 
a Schulz distribution 1311 with q — 5 at a packing fraction 
r/ = 0.4. The symbols are MC simulations.— The lines are 

SPT (•••), BGHLL ( ), VS ( ), eCS2 ( ), and 

eCS3 (— ). 



FIG. 6: (Color online). Plot of the contact value gn(cri) 
as a function of the packing fraction 77 for the binary hard- 
sphere mixture x% = 0.005, 0-2/0-1 = 0.2. The symbols are 

MC simulations^ The lines are SPT (•••), BGHLL ( ), 

VS ( ), eCS2 ( ), and eCS3 (— ). 



fact that it is the only previous theory consistent with 
the two ways of computing the pressure of the system. 
It is clear that the best overall performance in the com- 
parison with the simulation data is given by <?£, CS3 , not 
only qualitatively but also quantitatively. Interestingly 
enough, <7^ CS2 and g^ s also do a good job in the cases of 
mixtures 1 (Fig. QJ and 4 (Fig. Q}. In general, one can 
see that the SPT overestimates the contact values, except 
for high values of a, and that the BGHLL prescription 
underestimates them. Also, although not included in the 
figure to avoid overcrowding, there is a very poor agree- 
ment of both <?^ CS1 and <?^ Y (which are linear in z) with 
the simulation data, both being underestimations. 

While the main interest of our present formulation was 
geared towards the polydisperse system near a hard wall, 
Eq. 125fl . it should be clear that our new proposal for the 
contact values of the rdf also applies to the bulk fluid, 
Eq. (|24|l . In particular, if the diameter distribution is dis- 
crete, the replacements a — > Cj and a' — > <Jj in Eq. i|24|) 
yield the contact values gij(crij), where i and j denote the 
i-th and j-th species, respectively, and oy, = (<7j + <7j) /2, 



with (Tfc denoting the diameter of a sphere of species k. 
Just to illustrate the kind of results our proposal pro- 
duces, we consider two examples of binary hard-sphere 
mixtures. In the first case, we show in Fig. a plot of 
the different theoretical predictions of gii(ci) as a func- 
tion of the packing fraction r), for a mixture having a 
mole fraction of the large spheres x\ = 0.005, and a size 
ratio (T2/cri = 0.2, together with the simulation results 
of Cao et a/iSi On the other hand, in Fig. [7| and for a 
binary mixture with oil&\ =0.1 and two values of the 
packing fraction (77 = 0.3 and 77 = 0.4), we display the 
behavior of g\\{a\) as a function of the mole fraction of 
species 1 derived from the various theories and the results 
of simulation by Lue and Woodcock 22 and Henderson et 
al. 2 - It is again clear from Fig. 0that gi\ SZ {ci) gives the 
best performance for this mixture. As far as the behavior 
with respect to the mole fraction of the large spheres is 
concerned, as already noted by Henderson et al.?2, the 
approximation of Viduna and Smith— does a very good 
job, especially for 77 = 0.4, but it goes wrong for small 
values of x\. Figure {7\ also indicates that our proposal is 
certainly the best for 77 = 0.3 and that it also accounts 
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FIG. 7: (Color online). Plot of the contact value <?ii(crii) as a 
function of composition of the large spheres xi for two binary 
hard-sphere mixtures with 02/01 =0.1 and packing fractions 
77 = 0.3 (top panel) and r\ = 0.4 (bottom panel), respectively. 
The symbols are molecular dynamics (MD) simulations!— 

The lines are SPT (•■■), BGHLL ( ), VS ( ), eCS2 

( ), and eCS3 (— ). 



tion, g(a,cr'), and of the wall-particle correlation func- 
tion, g w (a), of a hard-sphere fluid mixture with an arbi- 
trary size distribution. The proposal relies on a kind of 
universality assumption by which, once the packing frac- 
tion is fixed, for all pairs of like and unlike spheres the 
dependence of the contact values on the diameters and 
on the composition is only through a single dimensionless 
parameter and holds for an arbitrary number of compo- 
nents. It also makes use of the point particle and the 
equal size limits and of the internal consistency between 
the usual virial route and the hard-wall limit route to 
derive the pressure of the system. As a consequence, the 
contact value g p of the rdf of a one-component fluid is 
required as the only input, thus making the formulation 
to be both simple and rather flexible. 

The merits of this proposal have been assessed by com- 
paring the contact values themselves and the correspond- 
ing compressibility factors with other theoretical devel- 
opments and with recent MC simulation results both for 
polydisperse hard-sphere fluids at a hard wall and for bi- 
nary hard-sphere mixtures (discrete size distribution) . It 
is fair to say that the new proposal with the CS expres- 
sion for g p gives the best overall performance. Also it is 
clear that (i) two different approximations g(cr, a') for the 
contact values can yield the same compressibility factor 
Z and (ii) a fortunate cancelation of errors can make a 
poor approximation for g(a, a') to lead to a reliable Z. 
Examples of the first effect are provided by the approx- 
imations PY and ePY3, which differ at the level of the 
contact values but share the same EOS, and, similarly, 
by the approximations BGHLL and eCS3. An example 
of the second effect is represented by the eCSl approxi- 
mation, which yields a very accurate EOS^ even though 
the associated contact values are only qualitatively cor- 
rect. 

In previous work of ours we have attempted to provide 
expressions for the contact values of the rdf and the com- 
pressibility factors that are valid for any dimensionality 
dM^& The present proposal does not fulfill such a condi- 
tion. It will also work for d = 1 and d = 2, but it is not 
directly generalizable to arbitrary dm* 

Apart from the EOS, contact values of the rdf may 
also be useful in other contexts. For instance, they are 
required as input in the rational function approxima- 
tion approach to the structural properties of hard-sphere 
mixtures. 25 We plan to use them in connection with this 
problem in the near future. 



correctly for the sharp rise observed at the smallest val- 
ues of the mole fraction of species 1 for both packing 
fractions. 

IV. CONCLUDING REMARKS 

In this paper we have provided a new proposal for the 
contact values of the particle-particle correlation func- 
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